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Abstract. We study the thermal and electric transport of a fluid of interacting Dirac fermions as they arise in single-layer 
graphene. We include Coulomb interactions, a dilute density of charged impurities and the presence of a magnetic field to 
describe both the static and the low frequency response as a function of temperature T and chemical potential /i. In the critical 
regime fl <T where both bands above and below the Dirac point contribute to transport we find pronounced deviations from 
Fermi liquid behavior, universal, collision-dominated values for transport coefficients and a cyclotron resonance of collective 
nature. In the collision-dominated high temperature regime the linear thermoelectric transport coefficients are shown to obey 
the constraints of relativistic magnet ohydrodynamics which we derive microscopically from Boltzmann theory. The latter also 
allows us to describe the crossover to disorder-dominated Fermi liquid behavior at large doping and low temperatures, as well 
as the crossover to the ballistic regime at high fields. 
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INTRODUCTION 

The quasiparticles in single layered graphene obey the 
massless Dirac equation [1, 2, 3]. At finite temperature 
and moderate doping, they form an interacting electron- 
hole plasma of ultra-relativistic fermions whose trans- 
port properties differ significantly from a conventional 
Fermi liquid. Indeed, it has been argued that intrinsic 
graphene at T = realizes a system at a quantum critical 
point which controls the transport properties in a window 
of temperature and chemical potential given by T > 
[4, 5]. This criticality is reflected in the inelastic scatter- 
ing rate being proportional to a^T (we set kg = h = 1 
throughout), /. e. , solely depending upon temperature. 
Here, a is the fine structure constant characterizing the 
strength of Coulomb interactions, a = e^/ ffv, where K 
is the dielectric constant of the adjacent medium and v is 
the Fermi velocity of the linearly dispersing quasiparti- 
cles. 

In this paper we show on the basis of a microscopic 
Boltzmann equation approach that at low frequencies 
and high enough temperatures graphene behaves like a 
plasma of charged particles and antiparticles obeying 
relativistic magnetohydrodynamics [6]. The temperature 
has to be high enough such that the inelastic scattering 
rate T^J is bigger than the driving a.c. frequency CO, and 
the scattering rates due to impurities, T^^, or weak mag- 
netic fields, Tg ' — v^eB/cT. AH the frequency depen- 
dent thermo-electric conductivities turn out to be com- 
pletely fixed by thermodynamic quantities, the covari- 
ant equations of motion and a single transport coefficient 
Gq{ij./T) which assumes a universal, entirely collision- 
dominated value in the critical regime /x/T ^ 1 of a 



clean system [5, 7]. 

A logarithmic divergence in the collinear scattering 
amplitude, occurring generally in two-dimensional sys- 
tems, allows us to obtain the transport coefficients to 
leading order in an expansion in the logarithm of the 
relevant infra-red cutoff, which in graphene turns out 
to be l/log(a). This approach also allows for an ele- 
gant description of the crossover from the critical regime 
ji < T to the Fermi liquid regime jj. ^ T, establishing 
the connection with existing results [8, 9]. To leading 
logarithmic order the collinear divergence results in an 
effective equilibration among quasiparticles moving in 
the same direction. This reduces the complexity of solv- 
ing the Boltzmann integral equation for the distribution 
function to the determination of at most three numbers 
that characterize the angle-dependent effective equilibria 
of collinearly moving particles. 

Very similar transport characteristics as in the criti- 
cal region of graphene is expected in other relativistic 
systems as they often emerge close to quantum criti- 
cal points with a dynamical exponent z = 1, such as 
the superfluid-insulator transition [10, 11] in a Bose- 
Hubbard model. Boltzmann transport theory has first 
been applied to such critical systems by Damle and 
Sachdev [12, 13], with the important difference that there 
only short-range interactions between the linearly dis- 
persing bosonic excitations were considered. Bhaseen et 
al. [10] have recently extended this approach, discussing 
interesting effects in the presence of a magnetic field. 
With minor modifications, the present general analysis 
applies to such systems, too. A detailed analysis of Boltz- 
mann magnetotransport of relativistic quasiparticles in a 
magnetic field was given in Ref. [15]. 
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RELATIVISTIC HYDRODYNAMICS 
FROM BOLTZMANN APPROACH 

If the Coulomb interactions are not too strong, a -C 1, 
a description of the system in terms of quasiparticles 
is possible. In graphene, this is certainly the case if the 
sample is adjacent to a dielectric with Jf :§> 1, but also 
holds generally at low enough temperatures, since the 
couphng a is marginally irrelevant and flows logarith- 
mically to zero as T ^ 0, see [16, 17]. This allows us to 
analyze transport within a quasiclassical Boltzmann ap- 
proach, solving for the distribution function fk^x^,t) in 
the presence of driving fields which vary slowly in space 
and time. 

The low energy excitations belong to two valleys 
around two inequivalent Dirac points K,K'. Within each 
vaUey the fermionic quasiparticles form two bands with 
linear dispersion e^ x = Xvk, where A = ±1 refers to the 
band of electrons and holes of the imdoped material, re- 
spectively. 

In equilibrium the distribution function is 



^0 



1 



exp[(ek,A-M)/7^] + l' 



(1) 



In the presence of an electrical field E(f), a temperature 
gradient VJ, and a perpendicular magnetic field B = B&^, 
the distribution function obeys the Boltzmann equation: 

d d 



E + v,,,x- 



= -^coii[A,r,k,?|{/}]. (2) 

Here ^cou[A,i",k,f |{/}] denotes the colhsion integral 
due to Coulomb interactions and impurity scattering, and 
y^ X = ^^v..X = ^k£k,A' denotes the quasiparticle velocity. 
Assuming small deviations from equilibrium the above 
can be rewritten as 

(dt + y (vk.A X e,) • ^ (r,f ) + Aoii[{/}] (3) 

= T [vk,;t •F^0^(k,A) +Vka •F^.?)^(k,A)] , 

where the RHS describes the driving terms due to small 
electric fields, and a slowly varying temperature, as char- 
acterized by the dimensionless vectors = and 
F^ = —vVT/T^. Further, we have defined 



<^.^(k,A) 



1 



(4) 



In linear response we can analyze each driving term 
separately. For the corresponding driving force we will 
simply write F below. In linear approximation we may 
parameterize the deviations from local equilibrium as 

/M(r,0-4',A(nr)) -4;,) X (5) 

xF- [yk,x8\\{^,k)+\xxez8±{^,k)] , 



where ^[| x are dimensionless functions of A and the 
modulus of k = vk/T only, while the dependence on 
time is understood to be impUcit. The electrical and 
the heat current associated with such a deviation from 
equiUbrium are given by 



Q/T J 2 V \ F(g|||0^)+e,xF(^x|</>O 

where the factor 4 accounts for both the spin and vaUey 
degeneracy, and the inner product 

{gl\g2)^Zj-0j2SliU)g2{U). (7) 

has been defined. 

The linearized Boltzmann equation can be cast in the 
generic form 



(6) 







where ^ = '^ + 1 



iO.. 



(8) 



(9) 



The Unear operators and ^ are the two components of 
the collision operator J^coU describing scattering due to 
inelastic Coulomb collisions and impurities, respectively. 
Explicit expressions have been given in Refs. [5, 15]. The 
dynamical term Q, acts Uke 

[i2^](k,A) = (Of^„il-jf,^)g{k,?i), (10) 

and arises due to the time derivative in (3). Finally, ^ 
describes the deflection due to the magnetic field, and 
consequently, it is the only term which mixes and 
gL- Typical eigenvalues of these operators correspond 
to the inelastic and elastic scattering rates t^":', t^' , the 
a.c. frequency CO, and the typical cyclotron frequency 
of thermal particles Tg^' = CO^ ~ v^eB/cT. One can 
show [18, 19] that the linear operators defined above are 
Hermitian with respect to the inner product (7). 

Here we concentrate on the regime of high temper- 
atures, low frequencies and weak fields such that 
dominates among the above rates. 

In the absence of a magnetic field one easily finds from 
the above the conductivity, J = OctE or the thermopower, 
Q = UxxE, with 



(y„ = 2j{(^''\^-'\(^^).. 



>• 



(11) 



(12) 



In the presence of a magnetic field this generalizes to 
J = CxiE + <7q,(E X Cz), wifli 



+ 



-1 d2^-l \aE 



ir), (13) 
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The above formalism is completely general and does 
not use any particular properties of graphene. In the 
sequel we will analyze a remarkable simplification which 
arises in the hydrodynamic regime ind = 2. Further, we 
will show how relativistic hydrodynamics emerges in the 
case of a massless linear dispersion relation. 



Divergent forward scattering in 2d 

It has long been noticed that interactions ind <2 lead 
to singular forward scattering between particles with a 
linearized spectrum [20, 21]. The physical reason for a 
diverging cross-section is that collinearly moving par- 
ticles share the same velocity to the extent that non- 
linear corrections to the dispersion can be neglected. 
Thus the interaction time in forward scattering processes 
diverges, leading to a singular amplitude. While in higher 
dimensions this effect is irrelevant, being suppressed 
by its small phase space volume, it emerges as a log- 
arithmic singularity in c/ = 2. The singularity is usu- 
ally cut off by non-linearities in the spectrum, which in 
graphene are predominantly due to renormalization ef- 
fects of Coulomb origin. In the case of RPA-screened 
Coulomb interactions in graphene there is another cut- 
off mechanism due to the vanishing of the full Coulomb 
propagator in exactly collinear scattering processes. Both 
mechanisms cut off the logarithmic divergence at small 
scales and lead to an extra factor of log(l/a) in the 
typical inelastic relaxation rates, i.e., {g\'^\g) I {g\g) ~ 
a2log(l/a)r. 

The strong forward scattering has an interesting phys- 
ical consequence: It leads to an effective equilibration 
among collinearly moving particles. In other words the 
distribution function will tend to a Gibbs distribution 
when it is restricted to quasiparticles with a particular 
direction of motion, yy^ x = e^, where is a unit vec- 
tor in the plane. However, the parameters of these one- 
dimensional quasi-equilibria will in general depend on 
the angle of motion with respect to the appUed field. As 
always, the equilibrium parameters (T, are conju- 
gate to the conserved quantities, i.e., energy, momen- 
tum and charge. However, for massless relativistic par- 
ticles, there is a further pseudo-conservation law which 
holds for weak interactions (a <C 1). In this limit, we can 
restrict ourselves to the Born approximation and eval- 
uate the inelastic collision term with Fermi's Golden 
rule. Within this approximation, the phase space for 
two-particle processes (1+2^3+4) in which the to- 
tal chirahty changes (Ai -h A2 ^ A3 + A4) vanishes ex- 
actly [13, 5]. In other words, to leading order in a^, 
the particle and hole number is separately conserved. 



which leads to an extra pseudo-equilibrium parameter, (p. 
We may call it a gravitational potential, since it couples 
equally to electrons and holes with positive energy. ' 

For a vectorial driving term (an electric field or a 
thermal gradient) we thus expect 



1 



exp 



ek,A-^"[*k,A]-k-(j"+^M[*Ml)-^^'^'[*k,A] 



+ 1 

(15) 
5f" 



where we have taken into account that 5r, 5m, 
can depend on the driving force F only hke 5r[v] = 
F-v5f, etc.2 

Note that in the case of a Unear dispersion ;^ = 
Xkv, the perturbation 5r is redundant since it can be 
expressed as a combination of two perturbations with 
8u/v = —5n/^i = 8T /T. Below we will thus concen- 
trate on the three soft modes ^ 



g(i{X,k) = Xk, 

g2{x:k) = X. 



(16) 

(17) 
(18) 



To logarithmic accuracy we can restrict the deviations 
from equilibrium to the above ansatz, since these modes 
are the only ones with a relaxation rate of order a^T (or 
smaller), while all other modes decay with a rate of order 
^(a2log(l/a)r). 

Notice that the response functions (p^-^ (4) are simply 
linear combinations of go and ^1, multiplied by (1 — 

f^x)- basis of soft modes is thus perfectly suited to 
capture the essentials of thermo-electric response. 

In a general two-dimensional Fermi liquid with arbi- 
trary interactions, one may still linearize the spectrum 
around Ef. However, the non-linearity in the dispersion 
of typical thermal excitations will be of order T /Ep, 
which now cut off the divergence, leading to an extra 



' In this paper we only consider perturbations E and Vr which do not 
change the overall number of particles and holes, but only redistribute 
them between different wavevectors. We thus do not have to worry 
about the very long relaxation time ■t^h ~ 'Cee/«^ which is needed to 
equilibrate the relative number of particles and holes. The situation 
is different, however, for scalar perturbations which do alter the total 
chirality [14]. 

^ For a scalar driving force such as V ■ S relevant in the case of bulk 
viscosity, the prefactor F ■ v in (16) is absent, while in the case of a 
tensoriaJ perturbation such as dui/dxj, relevant for shear viscosity, the 
prefactor has to be chosen as Y,ijViVjdui/dxj. 

^ In the presence of a perpendicular magnetic field there are additional 
term analogous to the above, but with v — > ez x v, and corresponding 
coefficients dn^,S<j>^,Suj^. 
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factor log (iiF / r ) as noticed in Ref . [20] . However, those 
authors concluded incorrectly that this logarithm would 
also show up in the transport relaxation time. Instead, 
we see here that there are always the above soft modes 
which decay with a slower rate without logarithmic en- 
hancement, because they represent deviations preserving 
unidirectional equiUbrium which kill the logarithmic di- 
vergence in the scattering rate. 



Recovering relativistic Hydrodynamics 

An entirely different approach to transport in the 
collision-dominated regime is provided by hydrodynam- 
ics. In the present case the relativistic nature of the fluid 
has to be taken into account, and one should analyze the 
equations of motion 

dpJP = 0, (19) 
dpTP^' = F«yjy+T;^^SpoT^", (20) 

where the stress energy tensor and the current of a rela- 
tivistic fluid are given by 



(21) 
(22) 



Here is the three-velocity of the energy density. We 
have also included a phenomenological momentum re- 
laxation rate due to impurities, '^"j,, which expUcitly 
breaks the Lorentz invariance of course. The correctness 
of this procedure for weak disorder will be confirmed be- 
low. 

The Reynolds-tensor t''^ describing viscous forces 
can be neglected for the response in the long wavelength 
hmit. The contribution v'^ to the electrical current cor- 
rects for the presence of energy currents that do not carry 
charge. The form (23) can be inferred by generalizing 
arguments by Landau and Lifshitz [22] for B = 0, in- 
voking covariance and the second law of thermodynam- 
ics [11, 6], which requires that the divergence of the en- 
tropy current be always positive. The argument leading 
to (23) is valid at small fields B, where it is allowed to 
linearize in B. However, it leaves undetermined the phe- 
nomenological conductivity parameter Gq > 0. 

The hydrodynamic approach allows us to obtain all the 
low frequency thermoelectric response by linearizing the 
equations of motions around equilibrium and reading off 
the response functions following the recipes of Ref. [23]. 
Only thermodynamic data and the single coefficient Gq 
enter into those expressions, which are discussed in de- 
tail in Ref s. [11,6]. This remarkable property of relativis- 
tic hydrodynamics, which strongly constrains the form 



of the response, is fully confirmed by the Boltzmann ap- 
proach, as we discuss now. 

To establish the connection between hydrodynamics 
and the microscopic analysis, it is central to notice that 
the "momentum mode" is an exact zero mode of the 
collision operator 'tago = 0. This is so because de- 
scribes the equilibrium of a system which moves at speed 
5m in the direction of the driving force. Since equilibrium 
is preserved by the translation invariant Coulomb inter- 
actions the perturbation can only decay due to impu- 
rity scattering, magnetic fields or the a.c. driving. In the 
hydrodynamic regime the latter all live on much longer 
time scales than the relaxation time of any other modes. It 
turns out that in this regime it is essentiaUy the dynamics 
of the momentum mode which governs the magnetore- 
sponse. 

The Boltzmann formaUsm provides an explicit expres- 
sion for the transport coefficient in the form 



where, in the absence of impurities. 



Omom — 



n {go\^\8o) n -icoe + P 



(24) 



(25) 



is the conductivity contribution due to the momentum 
mode. Here, £ is the energy density and P is the pressure 
of the free electron gas. In general amom diverges in the 
d.c. limit o ^ 0. However, at particle-hole symmetry 
(jj. = 0) the coupling to the momentum mode vanishes, 
(^^ |go) = 0, and Oq is seen to equal the d.c. conductivity 
in the absence of impurities. The interesting fact that this 
is in fact finite will be discussed below. To logarithmic 
accuracy one finds the result [5] 



Gq{ix = 0,(0 = 0) 



2e' 



0.121 



(26) 



For general H^O one can obtain a logarithmically accu- 
rate approximation for (7^ and Gq by observing that G^x 
is the maximum of the functional 



Q[g]^^[2{^^\g)-{g\ 



(27) 



and extremizing it over the hnear space spanned by the 
three soft modes ^o,,?!) ^2, Eqs. (16-18) discussed above. 
In full generality, one finds that Gq is a scaling function 
of ji/T which decays as a power law for jti » J [15]. 

The response functions Gxx, OCxx and the thermal con- 
ductivity Kxx have all the same structure in the absence 
of a magnetic field. There is a contribution from the 
momentum mode, which generally diverges in the clean 
d.c. hmit, and in addition there is an extra piece propor- 
tional to Gq which is contributed by the Coulomb relax- 
ing modes. An explicit calculation [15] shows that the 
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relative deviations between hydrodynamics and the mi- 
croscopic Bohzmann calculation are of higher order in 
^(oj Tee, Tee/ limp), which are indeed small in the hydro- 
dynamic regime. The same statement holds also for the 
case with magnetic field, where relative corrections are 
suppressed by factors of is)- 

Impurity scattering 

The momentum relaxation rate which was added in 
the hydrodynamic formulation (20) in an ad hoc manner, 
turns out to reproduce the leading disorder effects of 
the Boltzmann theory provided that the relaxation rate 
is taken to be the inverse of the elastic lifetime of the 
momentum mode, 

1 ^ (gol^lgo) ^ "imp 

(go|ii/»b) )c2max[r,|M|]- ^ ^ 

The final estimate applies to charged, unscreened impu- 
rities of density wimp- The lowest order corrections to the 
clean hydrodynamic results are simply obtained by re- 
placing CO — * OJ -h JTjjjjp. In particular the contribution of 
the momentum mode to the d.c. conductivity takes the 
generalized Drude form 

^^mom ' ~I~ Timp ; 7» ' (^9) 

h ^ e + P 

where the band mass m is replaced by (e + P)/{pv^). 
Crossover in the electrical conductivity 

The above described approximation of restricting to 
the soft modes lends itself naturally to the description of 
an experimentally relevant situation [24]: the crossover 
observed when the system is driven from the critical, rel- 
ativistic regime (/i <C T) to the conventional Fermi liq- 
uid regime (jj. ^ T) by adjusting the doping level, see 
Refs. [4, 15]. In the Fermi Uquid regime we can make 
connection with existing results [8, 9]. It is worth notic- 
ing that the asymptotics for jj. <^ T and jU 3> T is al- 
ready well described by only two modes, go and gi. In- 
deed, the "chiraUty mode" g2 either contributes very little 
to the current since it is hardly excited by the electrical 
field at <C r, or it is essentially identical to the "chem- 
ical potential mode" when ^^T, where essentially 
only one type of charge carriers with a definite sign of 
A is present. The crossover is thus expected to be rea- 
sonably well captured by this two-mode approximation, 
which was carried out in detail in Ref. [15] 

For low doping, in the critical regime jX <^T, the 
inelastic scattering time behaves as 

= a^T (30) 



reflecting the quantum criticaUty of the system [4, 5]. 
Upon doping, the Coulomb scattering rate decreases and 
ultimately behaves as xz} oc for ii^T if screen- 
ing is taken into account. This is the standard inelastic 
scattering rate in a Fermi liquid. The relaxation due to 
impurities behaves very differently as a function of T and 

as seen from Eq. (28). The increase of the scattering 
rate with decreasing temperature reflects the increased 
crosssection for Coulomb scattering at low energies. The 
crossover from a collision-dominated, critical conductiv- 
ity where a w (7q(jU = 0) ^ /ha? to an impurity lim- 
ited conductivity which is dominated by the contribution 
from the momentum mode (29) occurs when the density 
of thermally excited carriers, pth ^ jifivf-, (holes and 
electrons) equals the density of impurities, nimp. 

Of course it is possible to go beyond the logarithmic 
approximation based on two or three modes, if the log- 
arithmic divergence is properly regularized. In particu- 
lar, it was shown in Ref. [15] that in the weak disorder 
limit one can absorb the effect of aU modes other than 
gOjgi into one (or two in the case of strong magnetic 
fields) effective frequency-dependent matrix element(s) 
entering all the response functions. This is a more pre- 
cise version of the hydrodynamic statement that at low 
fields the whole response contains solely Og as non- 
thermodynamic parameter. 



Criticality at the Dirac point 

The fact that Oxx is finite in a clean, undoped system 
{[i = 0) is a pecuhar feature which is also known from 
neutral plasmas. It is possible because at particle-hole 
symmetry, an electric field does not excite a momentum 
flow, which would not decay under Coulomb collisions 
and thus yield an infinite conductivity. Instead, the elec- 
tron and hole currents decay within the finite relaxation 
time Tee due to their mutual friction. Since there is no di- 
mensionless parameter in the problem except for a, the 
universal colUsion-dominated scaUng (26), a ~ e^/ha}, 
can be predicted simply from dimensional considera- 
tions [25]. 

It is also interesting to study other thermoelectric 
response functions at the Dirac point and their crossover 
to to Fermi liquid behavior at large doping. For ex- 
ample one finds an anomalously large Mott ratio, 
—axx/{da/dfi), in the regime ju <i; as well as strong 
modifications of the Wiedemann-Frantz law, given that 
the thermal conductivity diverges in a clean system as 

^ 0. In both cases the deviation from Fermi liquid 
behavior at low doping, is diminishing at T ~ /x, re- 
covering the usual Fermi liquid characteristics at high 
doping levels. 

EquaUy unusual behavior of clean undoped graphene 
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can be expected with respect to the shear viscosity, a 
quantity which is of great interest recently in the con- 
text of another relativistic fluid: the quark gluon plasma 
where the shear viscosity is found to be anomalously 
small, approaching a lower bound conjectured from 
a general holographic correspondence with black hole 
problems and their area law. The large critical scaling 
of Tgg' '--^ r in graphene also suggests an exceptionally 
small value of the shear viscosity at the Dirac point (since 
77 ^ Tee, as opposed to the Fermi liquid regime jj.^ T. 
This is indeed confirmed by expUcit calculations. [14] 



Magnetotransport and cyclotron resonance 

An interesting prediction of relativistic hydrodynam- 
ics is the appearance of a collective cyclotron reso- 
nance [6], whose existence is confirmed by the Boltz- 
mann approach. As we have mentioned before, the dy- 
namics of the momentum mode dominates most of the 
thermo-electric response, since it has the longest life 
time. In the presence of magnetic fields, impurities and 
a finite a.c. frequency, go is not a zero mode of (8) any 
more, but there is still a unique pair of modes correspond- 
ing to the smallest eigenvalue of that operator (the orig- 
inal pair of zero modes g|| = gj_ = go is split into two 
modes). 

Within first order perturbation theory, projecting the 
operators in Eq. (8) to the momentum mode, one expects 
a resonance in the clean material at 



(Oc = ± 



{goW\go) 



= ± 



eBpv^ 



{gom/co)\go) c{e+P)- 



(31) 



Notice the similarity to the expression = eB /cm for 
a standard cyclotron resonance. Again the band mass is 
replaced by (e + P) /pv^, and ©c is proportional to 
the density, which reflects the coUective nature of this 
hydrodynamic mode. 

Including weak impurities to first order perturbation 
theory, and the coupling of the momentum mode to the 
faster relaxing modes in second order in B, one finds 
an imaginary part to the low eigenvalue. The frequency 
pair for which the operator on the left hand side of (8) is 
singular, becomes 



(Opole«±«c-i(7+ Vp). 

where the coUision induced damping is given by 



(32) 



(33) 



The pole OJpoie shows up in all frequency-dependent re- 
sponse functions as a resonance ~ (co — cOpoie)^^- The 
fact that 7 has such a simple expression in terms of the 



coefficient (Jq defined above in (24), is a special fea- 
ture of the massless Dirac spectrum which ensures that 
^go 0= (^^. At = the above relation (33) is indeed 
easy to check in second order perturbation theory, using 
that (7g = f7 is given by (11). A general proof of (33) 
was given in Ref. [15]. It is also shown there that the 
property ^go °^ 0^ is ultimately responsible for the fact 
that the relativistic hydrodynamics emerges in the micro- 
scopic theory. 

The interaction induced broadening 7 is already pre- 
dicted by hydrodynamics [6], and is exactly reproduced 
by the Boltzmann formalism. However, the latter allows 
one also to go beyond small magnetic fields, an assump- 
tion that had to be made when deriving the relation 
(23). While there is still a collective cyclotron resonance 
due to the damped dynamics of the original momentum 
mode, the pole must be found from the solution of 



Det 



(Wpole) 



lole 



0. 



(34) 



In perturbation theory one can estabHsh that % [cOpoie] ex- 
ceeds the lowest order result (32), while (33) overesti- 
mates the damping 3 [o)poie]. A numerical solution of the 
above secular equation shows that at high fields, where 



^-1 



Tge^ this trend persists. [15] 



Interestingly, this behavior is qualitatively similar as 
in nearly critical, relativistic conformal field gauge the- 
ories which are exactly solvable thanks to the AdS-CFT 
correspondence. [26] 

Clearly, the magnetic fields need to be such that the 
Landau level spacing Ell = hv^2eB/hc is still small 
compared to the temperature. Otherwise, interference ef- 
fects and the Quantum Hall effect start to set in, which 
cannot be described on the level of a semiclassical Boltz- 
mann theory. However, at small a, this still leaves a large 
window a^T <Tg^ <T for the Boltzmann theory to ap- 
ply in a regime where the fields have to be considered 
large and the straightforward hydrodynamic description 
with the constitutive relation (23) breaks down. 

In conclusion, we have discussed the transport prop- 
erties of a weakly coupled fluid of massless Dirac par- 
ticles. The critical regime of low doping and high tem- 
peratures is predicted to exhibit several unusual char- 
acteristics, such as a universal collision-dominated con- 
ductivity, an anomalously low viscosity and a collective, 
interaction-damped cyclotron mode. We hope that these 
effects will soon be within experimental reach in very 
clean, suspended graphene as studied in Refs. [27, 28]. 

The fact that in two dimensions forward scattering is 
anomalously strong allows one to solve the Boltzmann 
equation in a simple, physically transparent way which is 
asymptotically exact for weakly coupled graphene, and is 
well adapted to describe the crossover to the Fermi liquid 
regime. 
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